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Abstract 

We derive the stochastic master equations which describe the evolution of open quan- 
tum systems in contact with a heat bath and undergoing indirect measurements. These 
equations are obtained as a limit of a quantum repeated measurement model where we 
consider a small system in contact with an infinite chain at positive temperature. At zero 
temperature it is well-known that one obtains stochastic differential equations of jump- 
diffusion type. At strictly positive temperature, we show that only pure diffusion type 
equations are relevant. 

1 Introduction 

The theory of Open Quantum Systems aims to study the time evolution of a small system T~Lq 
interacting with an environment £ , cf [T)l2l ll7l ll8|. Starting from an Hamiltonian description 
of the coupled system [2J [TTJ Q2], the evolution of the reduced system Hq is obtained by 
tracing over the degree of freedom of the environment. 

In the Markovian approach of open systems, the time evolution of the state of the reduced 
system is characterized by a semigroup of completely positive maps, with a typical generator 
called Lindblad generator, which gives rise to an ordinary differential equation called master 
equation [2 El [20]. 

In this framework, an active line of research, motivated by recent experimental applications 
in quantum optics and quantum communications, is focused on the description of quantum 
measurement ^ [Zl |B1 |S1 HDl US H21 US [IZl HS1 HD1 HS|- Basically, in order to avoid Zeno 
effect [17] . the measurement is performed on the environment. According to the postulates of 
quantum mechanics, this involves a random perturbation of the evolution of the state of Hq. 
The dynamics of %q is then described by classical stochastic differential equations, which are 
perturbations of the master equation in terms of white noise [61 [151 El El El QUI [HI 111 EE 
[29IM1E5]- Usually, these equations are called Stochastic Schrodinger Equations or Stochastic 
Master Equations and their solutions are called Quantum Trajectories (the name "stochastic 
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Schrddinger equation" is usually reserved for the evolution of the state of Hq in terms of pure 
states whereas stochastic master equations concerns evolution of density matrices). 

In the literature, most of the results concern models where the environment is at zero 
temperature and a lot of investments are still in progress to get right descriptions for model 
at positive temperature (even without measurement). Typically, at zero temperature, for the 
random evolution of density matrices we get two typical stochastic master equations of the 
following form. 



pt)dW t , 



1. A diffusive equation 

dpt = C{ Pt ) dt + (Cp t + Pt C* - Tr [p t (C + C* 

where W% is a one dimensional Brownian motion. 

2. A jump equation 

d Pt = c( Pt ) dt + ( ^pfc*] " pt ) " Tr [c pt C ^ dt ) ' (2) 

where (Nt) is a counting process with stochastic intensity t — > Jq Tr [C p s C*]ds. 

In the above expressions C is the Lindblad operator. Note that we can recognize from the 
form of these equations that they all are simulations of the master equation, for they are 
stochastic differential equations valued in the set of states of %q and on average they satisfy 
the master equation. 

More complicated models use jump-diffusion stochastic differential equations which are 
mixing of equations CD and © 0011 [22] (see Section [2J. 

Mathematically, there are three usual ways to justify these equations. A first approach 
is based on Instrumental Operator Processes connected with the notion of Operator Valued 
Measures and Quantum Markov Semi-groups analogy [21 [6]. The second one is based on 
classical stochastic differential equation theory and on the concept of a posteriori state [6l 
El El [TUl I2T] . The third approach often called Quantum Filtering is based on the formalism 
of Quantum Stochastic Calculus and the notion of input and output field in quantum optics 
[21 El U21 [131 HH H3 [16]. In this last setup, the evolution of the small system and the 
environment is modeled by Quantum Langevin Equations (also called Quantum Stochastic 
Differential Equations or Hudson-Parthasarathy Equations) [HI2J122]- These equations are 
namely driven by Quantum Noises [H [2]. Next, by adapting the classical framework of 
stochastic filtering, one can obtain appropriate stochastic differential equations driven by 
classical noises, which take into account the "incomplete" information of indirect observations. 

In [231 [2H [25] , an alternative discrete way has been developed. The approach is based on 
the model of Quantum Repeated Interactions which provides a "useful" discrete approximation 
model of quantum Langevin equations [31 H]- The setup of quantum repeated interactions is 
the one of the interaction of a small system T-Lq with an environment represented by an infinite 
chain (S)fcLo^fc- Moreover, the pieces of the chain are identical and independent quantum 
system, that is T~Lk = H for all k. Each copy T~L interacts with T~Lo, one after the other, 
during a time r. In this framework, an appropriate language of discrete quantum noise and 
discrete quantum stochastic differential equations is constructed. Next, it is shown that the 
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time continuous limit (r — > 0) gives rise to continuous models of quantum Langevin equations 
and continuous quantum noises. In this context, a discrete time model of indirect quantum 
measurements, called model of Quantum Repeated Measurements has been developed in |23t 
|2~H |2"5] . It consists in performing a measurement of an observable of after each interaction 
with %q. As in the continuous case, the measurement introduces a random perturbation of 
the evolution of %q, described by discrete stochastic master equations and discrete quantum 
trajectories. In the same spirit as [1], by considering the continuous time limit (r — > 0), the 
continuous models of stochastic master equations and quantum trajectories are recovered. 
Much beyond the approximation result, these discrete time models (without measurement [1] 
and with measurement [231 IS1 EH] ) provide a concrete and intuitive physical justification of 
the continuous models of quantum Langevin and stochastic master equations. 

It is important to notice that these four different approaches of stochastic master equations 
are crucially based on the "zero temperature" assumption. For example in the approach based 
on quantum stochastic calculus, one the main obstacle to consider positive temperature model 
is to describe the action of the heat bath on the small system and to derive adequate Langevin 
equations. 

Recently, in [5], the discrete approach of quantum repeated interactions [3] has been 
adapted to models with heat bath at positive temperature. This way, thermal quantum 
Langevin equations have been obtained and a clear justification of the action of the thermal 
bath has been presented. It is then natural to combine this approach with quantum repeated 
measurements in order to derive stochastic master equations for heat baths. 

In order to introduce temperature, we consider that each copy % is in a thermal Gibbs state 
at inverse temperature f3. The crucial point in all the different approaches at zero temperature 
is the fact the state of the environment is a pure state. This is clearly not the case for a Gibbs 
state. In order to get around this difficulty, we apply the G.N.S. representation of that state. 
This way, the Gibbs state of each copy H can be considered as a pure state in an enlarged 
Hilbert space. Hence, with this representation, the convergence result of [25] can be applied 
and stochastic master equations for heat bath are derived. Surprisingly, models of the form 
([2]) with counting processes disappears and only diffusive models remains. 

This article is structured as follows. 

In Section [21 we remind the discrete models of quantum repeated interactions and then 
quantum repeated measurements. Next, we recall the main result of [25j which gives the 
stochastic Master equations as continuous limits of these discrete models. 

In Section O we adapt the result of [25] for model with positive temperature. This is 
achieved by describing the G.N.S representation of the heat bath. Hence we obtain the 
complete description of stochastic master equations for a small system in contact with a heat 
bath and undergoing indirect quantum measurements. 

2 From Discrete to Continuous Quantum Trajectories at Zero 
Temperature 

In this section we recall the mathematical description of the quantum repeated measurement 
model, as developed in [25] . 





3 



2.1 Quantum Repeated Interactions 

The model of quantum repeated interactions consists in studying the interaction of a finite 
dimensional quantum system T-Lq in contact with an infinite chain T<£ = Hk, where 
n k = H = C N+1 for all k. 

We first need to make precise the definition of the countable tensor product T3>. To 
this end, let us define an explicit orthonormal basis of T<£. Let {Xq, X\, . . . , X^r} be an 
orthonormal basis of T~L = C^ +1 . Actually note that only the choice of Xq is relevant in our 
construction: it represents a reference state of % (for example a ground state). Put Xf to be 
the copy of the basis vector Xi but acting on the n-th copy of T~L. The orthonormal basis of 
T<J> is then made of those tensor products 

X I ® . . . ® X? ® . . . 

such that all the i n 's, but a finite number, are all equal to 0. 

On T~L consider the basic operators Oj,i,j = 0, . . . , N defined by a l ^X k = b^Xj (in Dirac 
notation a*- = \Xj)(Xi\). We dilate them as operators a^(k) on T<I> by asking them to act as 
Oj on the k-th copy of T~L and as the identity operator on the rest of the chain. 

The model of repeated interactions [1] is now described as follows. Each copy of H is 
supposed to interact, one after the other, with %q during a time duration r. Each elementary 
interaction between %$ and T~L is described by a total Hamiltonian 

H tot = H ®I + I®H R + XHj . (3) 

The operator Hq corresponds to the free Hamiltonian of the system T~Lo, the operator Hr is 
the free Hamiltonian of the system T~L, the operator Hj is the interaction Hamiltonian and A 
is the coupling constant. 

The basis {Xq, . . . , X n } is chosen to be the basis of eigenvectors of Hr, that is, with our 
notations: 

N 

H R = J2lia°i4. (4) 

i=i 

The interaction Hamiltonian Hj is chosen to be of so-called "dipole-type" : 

N 

Hi = Y^ {Ci ® a° + Ct <8) 4) . (5) 

i=l 

After a time duration r of interaction, the evolution of Hq <8> H is governed by the unitary 
operator 

U = e - irHtot . 

That is, in the Schrodinger picture, the evolution of states on T~Lq ® % is given by 

P ' — y U P U* 

and in the Heisenberg picture, the observables evolve as 

X ^ U*pU . 

Now, in order to describe the repeated interactions on the whole space T-Lq®T^, we consider, 
for each k £ N*, the unitary operator Uk which acts like the operator U on the tensor product 
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T~io®T~Lk an d like the identity operator on the rest of the space. For k being fixed, the operator 
Uk describes the fc-th interaction. The whole procedure is then described by the sequence of 
unitary operators defined by Vt = UkUt-i ■ ■ ■ U\. For example the evolution of an initial 
state p on Hq <g) T<!> after k interactions is given by 

p ' y V kP V£ . 



2.2 Quantum Repeated Measurements 

Now, we are in position to describe the model of quantum repeated measurements |23j. To 
this end, we need to specify the reference states of Ho and %. Let p denote the initial state 
of %q. For each copy of T~L, we consider the usual thermal Gibbs state at inverse temperature 

Pfi = Tr[e-e»*y (6) 

where Hr is defined in expression 0. In particular pp is diagonal, with diagonal elements 
that we shall denote by {/3q,/3i, . . . , f3jy}- 

In order to describe the indirect measurement of an observable A of 7i, we come back to 
the description of the first interaction in the space %q (£> %. After the interaction the new 
state of %q (8> H is \i = U{p (8> pp)U* . Now if A owns the spectral decomposition 

v 

A = Y J \Pi, 

the measurement of A gives a random result in the set of eigenvalues Aq, . . . , A p . In particular, 
the value Aj is obtained with probability 

P [ to observe A* ] = Tr [p I ® P % ] . (7) 

After having observed the eigenvalue A,, the state p is projected and becomes 

(I®P i )p(I®P i ) 



Tr[p (I® Pi, 



For i being fixed, the state p\{i) represents the new state of %q ®T~L after the first interaction 
and the first measurement. Usually, we are only interested in the reduced system T~Lo, that is, 
we shall consider only the partial trace 

pi{i) = Tt-h (pi(i)) ■ 

The state p\ is a random state which takes the values 

A(p) 

T± [d{ P )] 

with probability Tr\Ci(p)] respectively, where 

d(p) = Tr n [(I 8) Pi) U(p ® Pf} )U* (/ 8) P^] . 
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This way, the state p\ is the new reference state of T-Lq and we can consider a new interaction 
with a copy of T~L and a new measurement of A. By reducing on Hq, we then get a new random 
state p2 which satisfies similar properties as p\. Iterating this procedure we obtain a random 
sequence of state {p k ) on rlo- This sequence, called discrete quantum trajectory describes 
the random modifications of the state of T~Lq undergoing quantum repeated interactions and 
measurements. More precisely, we have the following proposition (cf [23 J. 

Proposition 1 The sequence (pk) is a Markov chain. More precisely, if pt = 0, the state 
Pk+i can take the values 

cm . = n 

TV [£,(#)]'* U '-'" P 
with probability pi{6) = Tr [A(0)], where d{6) = Tr w [I ® Pi U{6 ® p p )U* I® Pi]. 

In other words, the above proposition can be summarized as follows. Let 1^ denote the 
random variable which takes the value 1 if we observe the eigenvalue A, during the k-th. 
measurement and otherwise, we then have 

2.3 Quantum Trajectories at Zero Temperature 

The equation above is a discrete-time stochastic master equation. In the articles |23[ [2^1 12"5] 
the author computes explicitely the continuous-time limit of ([8|) at zero temperature. In the 
limit, he obtains the usual stochastic master equations describing continuous measurement 
experiments. In general, these equations are stochastic differential equations mixing diffusive 
and jump noises. We shall now recall these results. 

We now focus on the case where the temperature is zero, that is (3 = +oo. This way, we 
have ^ = for alH = 1, . . . , N and pp = a® = \Xq)(Xq\. 

If K is an operator on T~Lq (g) T~L and for the choice {^o 5 • • • > Xn} of an orthonormal basis 
for T~L, the operator K can be written as a (N + 1) x (iV + l)-block-matrix, with coefficients 

Kj = Tr n [(iQlXJiXiDK 

being operators on %q. With these notations, we have 

N 



t=0 

In particular the A(p) are easy to compute explicitely: let U = (U l k )o<k,i<N be the block- 
matrix representation of the unitary evolution U and let Pj = (pL)o<fc KN m t ne basis 
{X , . . . , Xi} (in block form we have I 8) p = (p kl I)o<k,l<N), we get 

N 

Ci{p) = Tt h [I® Pi U(p ® pp)U* I®Pi\=J2 Pkl Vk P (V?)* ■ ( 9 ) 

k,l=0 
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The convergence result of discrete quantum trajectories is based on the asymptotic as- 
sumptions described in [3|. These assumptions concern the unitary operator U. If we consider 
the time of interaction r being r = 1/n, the unitary operator U depends on the parameter 
n, that is, U = U(n) = (U l k (n))o< k j<N. In [3], it is shown that the operator process (VL t ]) 
satisfying 

V[ nt ) = U [nt] {n) . . . Ui(n) 

converges, non trivially, to a process (Vt), only if the coefficients Uj(n) obey certain normal- 
izations. The limit process (Vt) then satisfies a quantum Langevin equation describing the 
evolution of a small system coupled with a Fock space. 

Their asymptotic conditions concern the existence of operators U- such that for all (i,j) G 
{0, . . . , N} 2 we have 

lim n^(U)(n) - 5^1) = iA, (10) 



where e# = ^(5 Qi + 5 0j ). 

In the context of measurement, the expression 



implies that only the asymptotic of 



the terms U9(n) are relevant. In terms of total Hamiltonian ([3]), in [1] it is shown that these 
asymptotics for the f/j's can be obtained by considering interaction Hamiltonian Hj of type 
([5]) and by considering the coupling constant A = \/n. In that case, they obtain 



-iC k 



Now with ()10p . we are in position to express the main result of [25J which links discrete 
and continuous quantum trajectories. To this end, we introduce functions (when it has a 
meaning) defined on the set of states: 



9i(p) 



N 



iu=i 



Tr 



N 



E AiL° kP (L 



0\* 

1 1 



fcj=i 



Vi(p) 

hi(p) 



Tr 



N 



k,l=l 



N 



J2{pi L° k P + ph k p(L° k r)-Tr 



k=l 



N 



k=l 



N 



£(p) = L p + p(L r + Y J L° k p(L k r. 



k=l 



Theorem 2 Let A = X^f=o^*^ ^ e an observable ofH. As ^^Pi = I, without restriction, 
we can assume that Pq ^ 0. Let I = {i £ {1, . . . ,p}/Pqq = 0} and J = {1, . . . ,p} \ I. Let po 
be a state on %q and let (p n (t)) be the stochastic process defined from the discrete quantum 
trajectory (p k ) by p n (t) = p[ n t\- We then have the following convergence result. 
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• If J = 0, the process (p n (t)) converges in distribution to the solution of the stochastic 
differential equation 

Pt = Po + C(p s -)ds + y2 9i{Ps-)lo<x<vi( Ps -)[Ni(dx,ds) - dxds], (11) 
Jo i=1 Jo Jr 

where (iVj)i<j<jv are N independent Poisson processes on R 2 . 

• If J ^ $, the process (p n (t)) converges in distribution to the solution of the stochastic 
differential equation 

Pt = po+ j C{p s -)ds+ [ hi{p s J)dWi(s) 

+ / 9i(Ps-)lo<x<Vi{ Ps -) [Ni(dx,ds) - dxds] (12) 

i£l Jo JR 

where (Wj(i))o<i<jv are N+l independent Brownian motions independent of the Poisson 
processes (Ni)i<i< N . 

2.4 The 2-Dimensional Case 

In order to illustrate this theorem, we investigate the case where T-L = C 2 . In this situation, 
we get two different behaviours depending on the fact that p® = 1 or not. 

Indeed, in the = 1 we have J = and the case Pqq ^ 1 corresponds to J ^ 0. 

Furthermore the = 1 corresponds to a case where the observable A is diagonal in the 

basis {Xq,X\}, that is, of the form 

A = A ag + Ai a\ . 

The limit equation is then 

Pt=P0+ f t £(Ps-)ds+ f [ ( jSro'-^L -Ps-)* 



Tr[L?p s _(L?)* 

x 1 o< a: <'rr[L;p s _(i;)*] [N(dx,ds) - dxds] . (13) 

By putting L\ = C and by considering the process N t = f* f R l <x<Tt[L° lPs -(LO)*] N (dx, ds), 

we obtain the jump equation ([2]) mentioned in Introduction. Indeed, the process (Nt) is a 
counting process with stochastic intensity J * Tr[Lip s -(L^)*] ds. Actually, the expression (fT3"j) 
is a rigorous way to consider jump stochastic Schrodinger equations (see [21] ). 

The other case p® ^ 1 gives rise to a diffusive equation. For example, consider the case 
Pqq = 1/2 (the other situations are similar). The observable A has then to be of the form 

A=^(aUal + al + a\) + ^(al-al-al + a\). 
Hence, we get the limit equation 

Pt = Po + J*Co(p s )ds + J* (x?p s + p s (i:?)*-Tr[L?p s + p s ( J LS)*]p s )^dW 1 ( S ) 

+ jf* (L\p s + p s (L\Y - Tr[L\p s + Ps (L?rK) ^ dW 2 (s) . (14) 
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Note that by defining a Brownian motion W t = {s/2/2) Wi(t) - (V2/2) W 2 (t), we recover the 
diffusive equation (TTJ). 

In [23l[23], the equations ()13| 14j) are studied in details and the convergence from discrete 
to continuous trajectories is obtained (with different techniques than the more general result 

125]). 

3 From Discrete to Continuous Quantum Trajectories at Pos- 
itive Temperature 

All the results mentioned in previous section are based on the construction of [3] which makes 
heavy use of the fact that the reference state of % is a pure state. Indeed, this condition is 
strongly needed in order to define the countable tensor product ® n ^*'H and its continuous 
limit, the continuous tensor product <8>t e K+'%. 

When considering that the environment is made of a chain of systems H each of which in 
thermal equilibrium state 




we cannot directly apply their results. The idea here follows the one developed in [5], that 
is, we take the G.N.S. representation of the state pp. This way, the state pp becomes a pure 
state, but on a larger state space. 



3.1 The G.N.S Representation of the Heat Bath 

The G.N.S representation of (H,pp), also called cyclic representation, is described as follows. 
At positive temperature, since the state pp is faithfull, it defines a scalar product on %' = B(H) 
by 

(A, B) = Tr [pp A*B] , (15) 

for all (A, B) € W . 

For all A £ T~L' , we denote by tt(A) the linear map from %' to W defined by 

n(A)B = AB, 

for all B € %' . The linear map tt from %' into B(H') is a representation of (T~L,pp), the 
so-called "G.N.S. representation". In particular we have 

(I,K{A)I)=Tr[ppA], 

for all A £ T~L' . This means that transported by tt the action of the state pp on a observable 
A is the same that the one of the pure state |/)(/| on tt(A). 

In order to compute the matrix coefficients of the operator n(U), we need to specify an 
orthonormal basis of Hi '. The only restriction on this basis is that it has to contain the 
reference state |/)(/|. 

The Hilbert space W is a (N + l) 2 dimensional space. We denote by Xq the identity 
operator. Next for % = 1, . . . , N, we denote by X\ the diagonal matrix with diagonal elements 
{vf,..., u^} such that 
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for all i,j = 0,...,N. Such operators can be constructed by extending the vector (1, . . . , 1) 
into an orthonormal basis of C^ +1 for the scalar product 

N 
i=0 

In order to complete the basis, we define X % - for i ^ j G {0, . . . , N} by 



a) 



Thus, we have construct an orthonormal basis {Xj,i,j = 0, ...,N} of %' for the scalar 
product (fT5j) . 

In this basis an operator K on %$ ® ~H is transported by ir as an operator ir(K) = 
(Kki)o<i,j,k,i<N where the coefficients are operator on These coefficients are given by 

K% = Tr n [(/ ® PP ){I ® ® Xj)] . (16) 

3.2 Asymptotics of U in the G.N.S. Representation 

Recall that we have defined the basic unitary interaction U as 

U = exp f-iL(H Q J + I (g, H R + y/n Hj^j . 

We are in position to describe tt(U) in the asymptotic way. Here, the translation of condition 
(|10p is the existence of operators LV, such that 



(17) 



J ki 

Jim n^i (u%-5 {id)m) I 



(18) 



where e[jg 



1,6' 



00 
kl 



M 
e 00 



1/2 and the others are equal to zero. 



We need to check that the unitary operator (|17p provides good asymptotic. Actually in 
our context of indirect quantum measurement, according to Theorem [21 we only need the 
expression of L^?- 

Proposition 3 The expression of the final relevant limit operators L?9 is given by 

N N 

L™ = iH + Y,PiHl + ^(PoCtCi + PiQCt), 



t00 
^k0 



1=1 



i=0 



Proof: In block form we have 

/ Hq v^C? y/EC$ ... 

V^Ci Ho+^I 

f^C 2 H + l2 I ... 



(19) 



H 



tot 



nC* N \ 





\Vn~C N 
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and hence the operator U(n) can be shown to be of the form (cf [3]) 



I-liH -ifal 

-YnClCt + o^) 



1 1 /-i r 1 * 



V 



l 1 fi /~1-k 



I-UH -i^ N I 

-lic N c% + o(^)j 



Now we can compute the asymptotic form of f7?P(n). Keeping in mind that in the ap- 
propriate basis B, the partial trace of an operator is the operator obtained by summing the 
diagonal blocks, we get 

= Tr n [I®ppX°UX°] = Tr n [I®ppU] 

= I ~ \ ( iffo + E ft-*' + 5 E (AC?Ci + AQC*) J + o (I) 
V i=l i=o / ^ ' 



(20) 



Let us stress that we have used 5^j=o^ = 1 to get the expression (f^O]) . Now if A; = and 
I 7^ 0, we have 



rrOO 
U 0J 



In the same way if I = and A; 7^ 0, we get 



-i=Tr w [(I 0^) (I<8> a ) C7] 
/A 

-i^v / A)C/ + (-p 



rrOO 
u k,0 



(21) 



(22) 



Actually the terms (|20t ETC] I22[) are the only terms which remain when considering the limit 
by applying (fTHj) . Indeed the other terms are expressed as 



rrOO 
u kk 



ifJL— 1 — / 1 N 

- E ^4(7^ + 2^*) + A)^ Ho/ + 2 E ^ 

V i=i V i=i 



+ ( - 



and if k 7^ 0, / 7^ and A; 7^ / 



rrOO 



1 1 

n2~' 



n 



Hence, by (|18p . they do not contribute in the limit. 



□ 
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3.3 Quantum Trajectories at Positive Temperature 

Before stating the equivalent of Theorem [2] with positive temperature, we need to be clear 
on how observables are transformed by the G.N.S representation. In particular, we have to 
describe 7r(J (g) P) when P is a projector. By the rule f)16[) . the coefficients Pj, of ir(I (g) P) 
are given by 

Pg = Ti:«[J®p /9 (I®X, fc )*J®P(7®A5)] 
= Tr„[J®(p^(Xf)*PXj)] 

= Tr[^(Xf)*PXj] J. (23) 

Define pjj = TV [ Pj3 [Xf )* P Xjj] and P' = (pg). Equation ([23]) means 

tt(P) = I <g, P' . (24) 

We are then in a similar situation as for zero temperature, but one has to notice a very 
important fact: the first coefficient p[jo ^ s now always strictly positive. Indeed, we have 

p™ = Tr [p p (X$y P X °] = Tr [ Pf} P] = Tr [ Pfj P*P] = (P, P) > . (25) 

This simple remark has an important consequence: there will be no jump contribution 
in the stochastic master equation for a heat bath. 

Indeed, let us consider an observable A = Y%=o an< ^ 7r (^-) = X^=o ^ ^ ® • ^ n 
Theorem [21 the jump contribution is directly connected to the set I = {i £ {1, . . . ,p}/p 00 = 
0}. In positive temperature, with the notation P' m = (jffA'rn)) the analogue of set I is the set 
I' = {i € {1, . . . ,p}/Poo{i) = 0}. Hence, the property ([25]) implies that I' = 0, which implies 
that there is no jump contribution. 

Consider the following functions defined on the set of the states: 



-Tr 



N 

-00\* 



k=l 

P l° (i) O + P ° i(i) L™ P + p™d) P (Ll° y + P ™(i) p(L™y 

N 



(26) 

C(p) = L^p + pL^UJ2{ L l> + P( L ^y + L 0kP + P(4°kr)- (27) 

fc=i 

Theorem 4 Let A = X/?=o ^i^i be an observable. Let po be a state on T-Lq. Let be 

the discrete quantum, trajectory describing the quantum repeated measurements at positive 
temperature. Let (p n (t)) be the sequence of stochastic processes defined for all t and all n 
by Pn{t) = P[nt](t)- Then (p n (t)) converges in distribution, when n goes to infinity, to the 
solution of the stochastic differential equation 

Pt = p + C( Ps )ds + y2 hi(Ps) dWi(s) , (28) 
Jo i=0 Jo 

where (Wi(t)) are N + 1 independent Brownian motions. 
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Proof: With the expression (|24j) . with the fact that the state pp is a pure state in the G.N.S. 
representation, we can apply directly Theorem [21 with the particular restriction we have 
mentioned above. This gives easily Equation ()28|) . □ 

With the explicit expression of the coefficients Equation (|28p can be made more 
explicit: 



Pt = Po + 



r* ( 1 N / 

/ - i[H Q , Ps ] - - {^{CiC k p s + Ps CtC k - 2C k Ps CI 

Jg V k=0 

1 N \ 

- 2 E (Afe^Cfc ^ + ^ C k C k - 2 Ct Ps C k )) ds 

k=0 / 



E 



771=0 

N 



N 



k=l 



Y,^{p^ k (m)C kPs -p^{m)p s Cl-Tv p™(m)C k p s -p™(rn)p s CZ]p s 

dW m (s). (29) 



o \/Poo( m 



k=l 



In this equation the expression 

A' 



C( P ) = -i[H , P]-\Y. (Po(<%C k p + P CtC k - 2C k p Ct 



1 N 

-Y J (h{C k Ctp + pC k Cl-2CtpC k )) (30) 



fc=0 



k=0 

corresponds to the usual Lindblad operator describing the evolution of a small system in 
contact with a heat bath at positive temperature. 

3.4 The 2-Dimensional Case 

As in Section 12.41 we want now to specialize the equation (|29p when H = C 2 and when 
considering particular observables. 

The first case is when the observable is diagonal, that is A = Xoclq + Aid}. In this case 
we have to compute vr(ag) and 7r(a\). As we have a[] + a\ = I and tt(I) = I, we just have 

to compute vr(a[j). Since Pqo(0) = Pqi(^) ana - Poo(0) = P?o(0)i we have only three terms to 
determine: pgg(0), p§?(0) and p§g(0). 
We get 



p°°(0) = Tr[ P/3 a°]=/3 



P8?(0) = Ttfo,(*y)*aM 



C^*„0 xrOl 



1 



pffi(0) = It^CJfo 1 )*^ 



/ft 
l 



Tr^a?a°] = 



(31) 



As a consequence the equation (|29|) for a diagonal observable becomes 



Pt = Po + / £(p s )ds 
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which is just the master equation for a heat bath, with no noise contribution. Hence, at posi- 
tive temperature, the repeated measurements of the observable A gives rise to a deterministic 
limit behavior. Recall that at zero temperature, the limit behavior was described by a jump 
equation. 

At positive temperature, we had already seen that the limit behavior should not involve 
jump contribution, but here, in addition we have no randomness at all in the limit. 

Let us now consider the observable A = Ao/2(a(j + a? + ciq + a\) + Ai/2(a(j — a? — ag + a\). 
As we have l/2(a§ + a? + clq + a\) + l/2(a{] — — + a\) = I, we just need to compute 
7r(l/2(a0 + a? + ciq + a\)) • With the previous computations, we only have to determine 7r(a°) 
and 7t(oq). In order to simplify the notations put R = a\ and 5 = aj. We have 

pOO _ oOO _ eOO _ c-00 _ n 

-"oo — -"-10 — °oo — °oi — u 
?oo _ nr coo 



Rqi — v A), Sio — y Pi ■ (32) 
Hence by ([311 El, for P = l/2(og + a? + + a{), we get 



P88CO) = |, Pg(0) = ^, rfg(0) = ^. (33) 
The equation ([29]) becomes 

Pt = Po + j* £{p s ) ds - J* ^ (ft (Cip, - p s C? - Tr[C lPs - Ps Ct} Ps ) + 

+ fa (Clps - p s d - Tr[Ctp s - p s Ci] Ps )^ dW (s) 

+ Si % ~T (f° (° 1Ps ~ PsC * ~ Tl[ClPs ~ psC * l]Ps ) + 

+ fa (c$p s - psd - Tv[Ctp s - Ps d\pM dWx{s) . (34) 

By defining a new Brownian motion W t = y/2/2(W (t) - Wi(i)) and by putting C = -id, 
the above equation becomes 



Pt = Po + f* £{Ps) ds + j I fa (Cp s + p s d - Tv[Cp s + p s d] Ps ) + 



+ fa (C*p s + p s C - Tr[dp s - PsC]ps)J dW s . (35) 

The equation ([35]) is then the equivalent of diffusive equation ^ (see Introduction) for the 
model with positive temperature. 
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